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Section I 

INTRODUCTION 


The Apollo Soyuz Te&t Project (ASTP) currently schedules an Apollo 
launch at 1500 EST on 15 July 1975 from Kennedy Space Center. The launch 
vehicle is subject to certain constraints, and the deadline for a "go" or 
"no go" decision for July 15 (or any subsequent date created by mission 
postponement) is at the previous midnight, i.e., 15 hours earlier. 

For mission planning, the effect of persistence in the weather events 
which cause operational delays is very important, so this study will take 
persistence into account in its presentation of conditional probabilities 
of launch conditions. Markov theory will be applied, and the results will 
be compared with empirical probabilities obtained directly from the data. 

Other studies of persistence in meteorological events have been conducted 
by Feyerherm and Bark (ref. 1), Williams (ref. 2), Hopkins and Robillard 
(ref. 3), Weiss (ref. 4), Eriksson (ref. 5), Brelsford and Jones (ref. 6), 

Smith (ref. 7), Gabriel and Neumann (refs. 8, 9), and Caskey (ref. 10). 
Additional material on runs is available in Gabriel (ref . 11) , Walker and 
Duncan (ref. 12), von Mises (ref. 13), Feller (ref. 14), and Wilks (ref. 15). 

A rather complete treatment of hypothesis-testing in regard to Markov 

chains is found in Anderson and Goodman (ref . 16) . 

\ 

As is well known, thunderstorms are a major impediment to launch operations 
at Kennedy Space Center in the summer. Two recent station studies of thunder- 
storms are the statistical investigations by Falls (ref. 17) and Neumann (ref. 
18). The former paper reaches the conclusion that thunderstorm events at Cape 
Kennedy are well represented by a negative binomial distribution. The latter 
paper discusses methods for predicting thunderstorms at the Cape, and it 
presents prediction equations obtained by nonlinear regression. 

The present study is divided into two main sections. Section II defines an 
ASTP unfavorable day and reveals the distribution of such "restrictions" over 
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the 15-year data period which is available. Then it analyzes the utility of 
past weather conditions in advancing forward a single day with conditional 
probabilities. Section III applies Markov chain theory and in particular, 
the Chapman-Kolmogorov Equations, to the data, comparing the results with 
experimental counts. Both first and second order Markov processes are 
investigated, and a statistical method is presented which predicts launch 
conditions on a dichotomous basis for mission planning. 

Finally, in Section IV a summary of the work is made and the conclusions 
are presented. The references are given in Section V. 
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Section II 

ASTP WEATHER RESTRICTION 
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2.1 DEFINITION 

An ASTP weather restriction is defined as the occurrence of any of the 
following conditions: 

a. Precipitation 

b. Thunderstorms with a cloud ceiling 

c. Cumulus cloud ceiling < 4000 feet 

d. Wind speed > 25 knots at 30 feet 

The record of hourly observations of weather elements at Cape Canaveral 
AFS from 1957 to 1971 is complete, and additional information is available 
on the peak wind speed each hour and the duration, location, and other 
characteristics of each thunderstorm. To treat just the 1500 EST observations 
would Ignore a substantial number of cases when restrictive conditions occurred 
during the afternoon and/or the weather at 1500 EST was nearly restricted. 
Therefore, the writers have adopted a broader definition of "unfavorable con- 
ditions" such that one or more of the four constraints is simply required 
to be reported at least once in the hourly record from 1200 EST to 1700 EST» 
inclusively. This makes the concept of "unfavorable" almost synonymous with 
the term, "marginal". 

A link can be provided between this definition and the probability of 
occurrence of restricted conditions at a particular hour. From all July and 
August data for the period 1957-1971, the relative frequency of restricted 
conditions at 1300 EST, 1400 EST, or 1500 EST, given that a restriction has 
taken place at least once between 1200 and 1700 EST, is 0.333, 0.390, or 
0.434, respectively. The 95 percent confidence limits on each value are 
+ 0.049. 

2.2 TIME DISTRIBUTION OF ASTP UNFAVORABLE DAYS 

Five-day moving averages of the frequency of unfavorable days for an 
ASTP launch are calculated for each July and August day and shorn in Figure 
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2-1. From the broken, visually-fitted curve, a rising trend is apparent in 
July with a reversal early in August. There are a few irregularities, but 
only the dip in mid-August is significant at the one percent level by Student’s 
t-test. In view of these results, the data are considered to be sufficiently 
homogeneous to permit the inclusion of August weather reports. This produces 
a total of 930 days of observation for the 15-year period under study. A 
test for cycles in the occurrence of unfavorable days discloses that the 
data are free of such periodicities. 


5 - DAY MOVING AVERAGES 



DAY 


Figure 2-1. FREQUENCY OF DAYS WITH UNFAVORABLE CONDITIONS, KSC, 1957-1971 

2.3 ANALYSIS OP A ONE-DAY ADVANCE 

The following question will now be considered: If an advance of one day 

Into the future la to be made In a prediction scheme baaed upon the past 
weather with Ite per Bistent nature, how many daya of peat weather can be pro- 
fitably used? 
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2.3.1 Theor*ticol Aspects 

Table 2-1 presents empirical and calculated values of favorable conditions 
for ASTP launch the following day, given present and previous days weather 
behavior. In the table P(F^) is t * ie probability of the next day being favor- ' 

able for the launch, independent of present or previous days weather* This 

\ 

is considered a zero order probability. P(Fj|U o ) is the probability of the ; 

next day being favorable for launch, given that the weather today is unfavor- j 

i 

able for launch. This is considered a first order probability. Similarly, | 

P(F 1 |U o F_ 1 ) is the probability that the next day will be favorable for launch j 

given that today is unfavorable and the day before was favorable. This is con- 5 

sidered a second order probability. The results of Table 2-1 are also plotted 1 

in Figure 2-2. ! 


Table 2-1 . EMPIRICAL PROBABILITIES AND CALCULATED PROBABILITIES 
OF FAVORABLE CONDITIONS FOR ASTP LAUNCH 



EMPIRICAL VALUE 

CALCULATED VALUE 

CALCULATED FROM 

P^) 

0.694 + 0.026* 



P(U,) 


0.306 + 0.026* 

P(F) 


0.788 + 0.028 

0.787 + 0.022* 

p< F i IV 

P(F,|U 0 ) 

0.483 + 0.051* 

0.481 + 0.064 

p < F .l F o> 

P < F ll F 0 F -1> 

0.823 +0.028 

0.807 + 0.020* 

P( F llF 0 U.,) 

P(F l' F o U -l) 

0.714 + 0.070* 

0.654 + 0.103 

P < F ll F o F -l> 


0.493 + 0.069* 

0.492 + 0.072 

P <W-1> 

p ( p l l°o F -l) 

0.473 + 0.077 

0.472 + 0.074* 

Pf^lV.,) 


The eubeoript8 indioate the order of days for a favorable (F) 
or unfavorable (U) oaee . The asterisks indicate the valuee 
used in eubeequent oaloulatione. 
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The empirical probabilities are obtained by counting the number c.f days 
that fit the class description (r) divided by the total number of day., in f e 
data set (n). For example, 

p(Fi!u o p. x ) * . (i) 

Thus by counting the n’lmber of days r in the total data set where a favorable 
launcn day is followed by an unfavorable day and a preceding favorable day, 
and then dividing this by *ne total number of days in the data set, n, an 
approximation to p(F^|U Q F_^) is obtained wnich Improves as n •> «, In this 
limit r/n -*• p(F^|U Q F_^) . Confidence limits are then calculated on p using 
Bayes Theorem as shown in Pratt, Raiffa and Schaefer (ref 19). There is some 
prior distribution of p assumed before the present data are utilized. This is 
called the prior distribution given by f T (p). Since all that is known is that 
p has a value between 0 and 1 and any value is equally likely, one can let 


f x (P> - 1 


( 2 ) 


The data, D , 
r ,n 

the data D , 
r,n 


is then obtained. Then the posterior distribution of p, given 

[f_(p|9 )], car. be obtained from Bayes Theorem 

F r gti 


f(D Jp)f (p) 

-«-l D r.n> V i 

t »n 


(3) 


The probability of r successes in n trials for a given value of p is given 
by 


f(D r,J p) " P r < l -P) n " r 


(4) 


This gives for the posterior distribution of p, efier the data have been 
utilized. 


f-' 


J /(1-P) n “ r dp 


( 5 ) 
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This is the well known Beta Distribution. The most likely value is given by 

the maximum of f_(p|D ), which can be found as 
F ! r,n 


<*Mp|d 


r iH 


dp 


) 


0 


( 6 ) 


which gives the most likely value of p as 


P 



(7) 


Then the 95 percent confidence limits, p and p. , are given by 

U L 

r*\l 

. fp(p|D r JdP “ .95 (8) 

* r » n 
P L 

It can be shown that the Beta distribution for large n can be approximated by 

2 

a normal distribution with a mean, u(p) of r/n and a variance, a (p) of 

— (1 )/n. Thus the normal tables can be used to find the 95 percent confidence 

n n r 

limits p u and p^. For the normal approximation this is given by 

p - ± ii!65 o(p) (9) 

L /n 

As shown by Cohen, in reference 20, the sero order, P(F ), first order, 

o 

P(F.|F ), and second order, P(F,|P F , ) , probabilities can all be obtained 

10 1 O “1 

from any four given empirical probability values. For example. Table 2-1 presents 
a set of calculated values based upon the empirical values for P(F 0 ), P(FjjF o ), 
P(F 1 |F o F_ x ), and U^). 


Some immediately derivable expressions are listed below as equations 10 
through 14. 

P(U ) - 1 - P(F ) (10) 

o o 


For a first order process, 

p(u x |Fo> - i - f(f x |f o ) 


(11) 


! ! I f 



I 

I 






and 



P(F ) 

p<r ilV •?( uT p ( u i | r o ) 

o 

While for a second order process , 

P (^i I P o p -l> ’ 1 - P(F ll P o F -l> 

and 

P(U l |U o U -l } " 1 - p ( F il U 0 U -i> 


( 12 ) 


(13) 


(14) 


Similarly* it hat bean ahovn by Cohen (ref. 20) that 


P ( Ui | F o U _ x ) 


P(F o ) [l-P(F o |F x ))- P(F^) P(F e |F 1 )(l-P(F 1 |F o F -x )] 
P(F 0 ) - P(F_ X ) ?(F o |F_ 1 ) 


This can be rewritten as 


POJJFo U_ x ) 


P(F 1 ) PtUjF^ - P(F a ) P(F 1 lF o ) P(U 1 [F o F^) 


P(F X ) [1 - P(F 1 |F o )] 


P(^ 1 lF (> ) ~ P(F l lF Q ) P(U 1 |F p F^) 


P ( U 1 lFo ) 


(15) 


This is equivalent to 

„„ I, , , P(P l'V P(U l |P o F -l> 
p < F l l F o F -1 ' FOIjl**) 

By Interchanging F and U* one can also obtain 

r(u,|u„) F(? 1 |U Q ij 


p « 1 l"o F - i > 

and equation (14) can be used to obtain 




(16) 


(17) 


p < r i|U 0 U -l> " 1 ’ P<U J°o U -l> 


08) 
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The remaining calculated values of Table 2-1 are complements of values 

found by equations (15) and (17). Differences between corresponding pairs of 

numbers in this table are but a few percent and they are always within the 

confidence limits. To illustrate the determination of confidence limits , suppose 

that P(FjU o ) is to be computed from PtF^jF^ by equations (11) and (12). Since 

P(Fi|F o ) is a random variable, its 95-percent confidence limits pj are given by 

equation (9) where p = p(F^|F Q ). Using the symbols < > to denote the averaging 

process, the confidence limits of the calculated quantity, P(U. |f ) are found 

2 1 o 

by deriving the variance a [P(U^|F q )] from 

° 2 E p ( U i 1 f q ) ] = "H - p(f 1 |f q )] 2 > - <[i - p(f 1 |f q )]> 2 

= <p(f 1 !f o ) 2 > - <P(F 1 |F o )> 2 - o 2 [P(F |F q )J (19) 

Since the variance is unchanged by subtraction, the confidence limits of the 

calculated value P(Uj|F q ) is equal to the empirically-determined value found 

for p(F |F ) . 

1 o 

Therefore, from Table 2-1 and equation (11), 

P(U. |F ) = (1 - 0.788) +0.028 
I o — 

= 0.212 + 0.028 (20) 

From equation (12), P^ju ) « K P(UjjF ) (21) 

where 


K 


2SI1 

P(U) 


Following the same procedure as before, 


a 2 [P(F 1 |U o ) ] - < [K P(U 1 |F q )] 2 > - <(K P(U 1 |F (> )1> 2 
- K 2 o 2 [P (U x | F q ) ] 


( 22 ) 


Thus the confidence limit is changed by a factor of K. Now from 
equations (20) and ( 7 1) ar.d Table 2-1, 
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n. 


"■< "'■m iwi f i v 1 







( 23 ) 


p < F ii V ‘ oil <°- 212 ± °- 028) 

= 0.481 + 0.064 

This result is shown in Table 2-1. Now since pCFjl^) and P(FjjF o ) can be 
found experimentally, either one can be calculated from the other. In making 
the choice of an empirical value, it is reasonable to select the value which, 
together with its corresponding calculated value, has the narrowest confidence 
limits. Such selections are marked with asterisks in Table 2-1 and are also 
plotted in Figure 2-2. 

2.3.2 Analysis of the Order of the Markov Process 

Figure 2-2 shows the probability of a favorable day for given conditions 
during the previous days. These results were taken from Table 2-1. 

If the present case were described by a zero order Markov process, the 
conditional probabilities, P(F 1 |F q ) or PCFjt^), would be equal to P(F 1 >. 
However, as shown in Figure 2-2 they are significantly different from the value 
of P(F X ), so that it is extremely unlikely that a zero order Markov Process 
would describe the present case. It is also shown in Figure 2-2 that the 
probability of the next day being a favorable day, given that the present day 
is favorable, will also depend on the previous day; that is, P(F 1 |F q U^) is 
significantly different from P(F 1 |F q F^) . This means that if the present 
day is favorable, the probability of the next day being favorable would depend 
on whether the previous day was favorable or unfavorable. 


Thus the condition P(- | F- > can be considered a second order process. 

However, in Figure 2-2 note that U^) and Ptfju, F^) are very close 

to P(F|U). This means that P(* |U) can be considered a first order process. 

The results for the probability of an unfavorable day are equal to 1 minus the 
results for a favorable day and are also shown in Figure 2-2. 

2.4 RUMS OF PERSISTENT ASTP LAUNCH CONDITIONS 

In counting the number of days in a run, all sequences which originate before 
July 1 are excluded and those which extend into September are included, bot 
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Figure 2-2. PROBABILITY (P) OF FAVORABLE (F) OR UNFAVORABLE (U) CONDITIONS THE 
NF’»T DAY, GIVEN 1 DAYS OF RECORD. THE 95 PERCENT CONFIDENCE LIMITS 
fi : ,c SHOWN FOR 1 < 2 
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2.4.1 Runs of Unfavorable Days 

The longest run of consecutive unfavorable days after a favorable day 
is one week long. This is shown In Figure 2-3, where relative frequencies give 
the empirical probability of runs of n or more days, with n ranging up to 7 days. 
Note that p(n>l) is the probability of one or more unfavorable days after a 
favorable day, and its empirical value of 0.212 is equal to 1-P(F^|F q ) in 
Table 2-1. 


The theoretical probabilities of the same runs are calculable from the 
basic empirical values found by counting cases. Thus the probability of n or 
more unfavorable days after a favorable day is given by 

P(U (n) F ) * P(F .) P(U |F .) P n_1 (U |U ) (24) 

—l -l o -l 1 o 

The equations for a chain of events is discussed more fully in Section III. 
Normalizing yields 

P(U (n) F_ 2 
y n “ P(F_ X ) P(UJ 


F -l> 


P n-1 (ujiy 


(25) 


Curve A (Figure 2-3) shows the result when the probabilities are evaluated 
from the data by equation (24), and it is in fine agreement with the empirical 
Information, at least to 0.1 percent. The utility of this graphical relation- 
ship is to read the probability of unfavorable ASTP launch weather for n or 
more days, given the present day is favorable. 


Since Curve A is straight on semilog paper for n >_ 1, then 


,-k<n-n ) 
o 


(26) 


Evaluating k directly from Curve A yields k ■ 0.658. Then for a one-day advance. 


P(U |U ) - e _k(1) - 0.518 

y n Vo 


( 27 ) 
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EMPIRICAL VALUES OBTAINED BY COUNTING ARE ALSO SHOWN AS INDIVIDUAL POINTS 


Figure 2-3. CUMULATIVE PROBABILITY (P) IN PERCENT OF THE OCCURRENCE OF N OR 
MORE UNFAVORABLE DAYS, AFTER A FAVORABLE DAY (CURVE A), EQ. (24), 
AND OF N OR MORE FAVORABLE DAYS, AFTER AN UNFAVORABLE DAY (CURVE B), 
BASED UPON 124 REVERSALS IN EACH CASE OUT OF A TOTAL OF 826 DAYS 
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This outcome for p (Uj_|U ) agrees closely with the empirical value (0.517) 
obtained by counting cases. A similar linear relationship on runs was also 
examined by Langley in reference 21 for wet periods in Montreal. 


2.4.2 Runs of FavorobU Days 

Using the same counting rule, the longest run of favorable conditions after 
an unfavorable day is found to be 25 days in length (Figure 2-3). Again, the 
probabilities have been estimated by relative frequencies counted from the 
data. The theoretical probabilities of n or more favorable days after an 
unfavorable day is given by 

P(F (n) U_ x ) = POJ^) P(F o |u_ 1 ) P(F 1 |F o U_ l ) P n_2 (F 2 )F 1 F q ) (28) 

Normalizing gives 

P(F (n) U_ x ) 

y 'n “ POJ^) P(F q |U_ 1 ) “ 

Curve B (Figure 2-3) shows this result when the probabilities in equation (28) 
are evaluated from the data. Agreement with the plotted points is fairly good, 
although there are trends away from the theoretical curve when the probability 
falls below 0.02. Again, the probability of persistently favorable launch 
weather for n or more days, given the present day has unfavorable conditions, 
can be read from the graph. 

Estimating k directly from Curve B yields a value of 0.203. Then for a 
one-day progression, 


P< F l|F o U^) P n <F 2 1 F F q ) 


(29) 


n+1 


P(F 1 |F o F^) - e~ k(1) - 0. 


816 


(30) 


This outcome for P(F^|F F_^) Is slightly less than the empirical value of 
0.825 obtained by counting. 
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Section III 

CONDITIONAL PROBABILITIES FROM THE CHAPMAN- K0LM060 ROV EQUATIONS 


In this section, the applicability of the Markov chain theory will be 
stressed, and the Chapman-Kolnogorov Equations (Breiman, ref. 22) will be 
uced to obtain conditional probabilities for up to 4 days following a base 
day with its identifiable favorable or unfavorable condition and available 
past weather record. 

3.1 MARKOV CHAINS OF ORDER ONE AND ORDER TWO 

The law of multiplication in probability theory states that the probability 
of a string of events E^, ..., E t+n is the product of n factors, 

P(E t , E t+1 . ... E t+n ) - P(E t > F(E t+l |E t ) P(E t+2 |E t+l> E^ ... 


P(E t+n' E t+n-l E t+1’ V 


A Markov chain of first order also has n factors, but it retains but 
one previous event as "given" in each factor, thus, 


P <V E t + 1’ ••• W - P( V P(E t + llV P <WWl> 


Analogously, a Markov chain of second order retains two previous events 
as "given" in each factor, so that 


P <V >w ••• W - p < E t> p<E t + il E t> p(E t« |E t + i’ V ••• 

E < E tJW’ E t+n-2 ) 


Many of the papers referenced in Section I which treat persistence in 
precipitation and in dry periods find that a Markov model, especially the 
first order chain, is an acceptable device to describe the data. 
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3.2 THE CH APM AN- KOLMOGOROV EQUATIONS 

The probability of passing from the initial state A i at time n to state 

A at time n+1 can be written as P(A Ja ). The probability of passing 
j j i f n 

from the initial state A^ at time n to a new state A^ at time n+2 is given by 
the Chapman-Kolmogorov equation for a first order Markov process as 


p(A k ,n+2 ,n+2 ^ A j ,n+l^ P(A j,n+l^ A i,n ) 


( 34 ) 


Thus the probabilities have been summed over all possible intermediate 
states Aj at time n+1. This process can be continued. 


P(A H,n+3l A i,n ) * £ P(A £,n+3l Vn+2^ P(A k,n+2l A i,n ) 


( 35 ) 


etc. The same type of analysis is applicable to second order or zero order 
processes . 


This result can be used to predict unfavorable or favorable conditions 
as far ahead as desired. The results for favorable conditions can be obtained 
from the unfavorable relationships by the following equations, 

P(F i |U o ) - 1 - P(UjU o ) (36) 

P(F i |F o F_ x ) - 1 - P(U 1 |F q F_ x ) (37) 

P (F il F 0 U_ x ) - 1 - P(U 1 |F o U_ x ) (38) 

These results can be used to predict the behavior in advance, given the 
present known weather conditions. The method of calculation is illustrated 
in Figure 3-1, and the outcome for U^, i ■ 1, 2, 3, 4 for July and August 
at KSC is shown in Table 3-1 and Figures 3-1 and 3-2. Experimental counts 
are included in the table for comparison with the theoretical predictions, 
and agreement is to within 5 percent absolute value in all of the conditional 
probabilities. This Indicates that the Chapman -Kolmogorov Equations are indeed 
applicable to this problem. 
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Table 3-1. CONDITIONAL PROBABILITIES OF UNFAVORABLE CONDITIONS 
FOR ASTP LAUNCH 



pi'W 


P(u, 1 

F 0 U -l> 

P(U,|F 0 

F -i> 

ith 

E* 

c** 

E* 

c** 

E* 

c** 

day 

1 

0.517 

0.517 

0.272 

0.216 

0.175 

0.193 

2 

0.410 

0.454 

0.267 

0.249 

0.242 

0.255 

3 

0.331 

0.372 

0.319 

0.270 

0.274 

0.284 

4 

0.306 

0.343 

0.245 

0.278 

0.278 

0.297 


"•Empirical value obtained by counting and use of the complementary 

relationship . 

** Calculated value from the Chapman-Kolmogorov Equations. 

P (Ui | u 0 ) - probability of the i^ unfavorable day given that today is un- 
favorable. 

I | F c U_i) - probability of the i th unfavorable day given that today is 
favorable and yesterday unfavorable. 

P(Ui | - probability of the ith unfavorable day given that today is 

favorable and yesteiday favorable. 

Unfavorable conditions mean that one or more ot the four launch constraints 
(see page 2-1) occurred at least once between 1200 EST and 1700 EST. 

Examples illustrating use of the tables : 

1. This afternoon was unfavorable for ASTP launch. The empirical and calculated 
probability that tomorrow afternoon will be unfavorable is .517. The empirical 
probability that the third afternoon from today will be unfavorable is 0.331. 

2. This afternoon was favorable, but the previous afternoon was unfaviiable. 

The calculated probability that tomorrow afternoon will be unfavorable is 0.216. 
The empirical probability that the 4th afternoon from today will be unfavorable 
is 0.245. 

3. This afternoon was favorable and the previous afternoon was favorable also. 
The empirical probability that tomorrow afternoon will be unfavorable is 0.175. 
The calculated probability that the afternoon after tomorrow (2nd day in the 
future) will be unfavorable is 0.255. 
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Figure 3-2. 


CONDITIONAL PROBABILITIES OF MARGINAL ASTP LAUNCH CONDITIONS 
IN JULY AND AUGUST AT CAPE KENNEDY, BASED UPON THE CHAPMAN- 
KOLMOGOROV EQUATIONS (P) AND UPON EMPIRICAL COUNTS (p), FOR 
1 DAYS IN THE FUTURE FROM DAY ZERO. 
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Section IV 

EXAMPLES OF PROBABILITY CALCULATIONS 

The values of Table 2-1 and their complements for unfavorable cases are 

sufficient to compute the probabilities of runs with no previous conditions. 

For example, if P(U U U ) is desired, equation (32) is applied to get 

o 1 2 

P(U U-U-) - P(U ) P(u.|u ) 2 - (0.306) (0.517) 2 - 0.082 
Oil o i o 

If P^F^F^) is desired, equation (33) is applied to yield 

P(F F ,F,) » P(F ) P(F_ I F F ) 2 - (0.694) (0.823) 2 - 0.470 

Oil O l IO 

In case a probability specified at a particular hour is sought, the 
appropriate link fro.2 subsection 2.1 is inserted. Thus, to find the 1500 EST 
probability of a restriction, given a previous day's restriction at 1500 
hours, use 

P(U 1500,ll U 1500,0 ) “ P(U o^ U 1500,0 ) P(U 1500,l^ U o ) 

“ P < U ol U 1500,0> P < U l'V P(U 1500,ll U l> 

- (1.00) (0.517) (0.434) - 0.224 

The unconditional probability of encountering two consecutive restrictions 
at 1500 hours is 


P(U 1500,1 ^OO.O* 


P(U 1500, 0 IV P(U o> P(U JV P(U 1500,llV 
(0.306) (0.434) (0.517) (0.434) - 0.030 
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Section V 

SUMMARY AND CONCLUSIONS 


In :he Apollc-Soyu* l: «. »' rogram, the declelon to launch must be made 
at least 15 hours before .. .nch time. Therefore, mission planning can be 
aided by a knowledge ci statistical relationships between the occurrence of 
inclement (and favorable) ASTP weather conditions in past time and future time. 
Markov theory has been applied in this study to elucidate these relationships 
in terms of conditional probabilities for Kennedy Space Center. 

The first forecasting problem investigated was the length of record of 
past weather which is useful to a prediction. Based upon the historical 
sequence of hourly reports for July and August from 1957 to 1971, relative 
frequencies of marginal ASTP weather were gleaned from the data and expressed 
as four empirical conditional probabilities from which other conditional 
probabilities up to second order were derived. The outcome is contingent upon 
the nature of the preceding weather. Thus, if afternoon weather for the current 
day has been unfavorable, the previous afternoon's reports hive negligible fore- 
cast value. On the other hand, if the afternoon weather for the current Jay 
has been favorable, the previous day's reports are important to the prediction. 
These results signify that first order and second rder Harkov cnalns, 
respectively, are operative. 

The second forecasting problem studied was the matter of runs of favorable 
or unfavorable launch conditions. Such runs were found to persist as long as 
25 days and 7 days, respectively. In the case of unfavorable present weather, 
there is a conditional probability slightly greater than 0.50 that inclement 
conditions will persist another day. On the other hand, the probability of 
one or more favorable days, givsn unfavorable present wssther, i® only about 
0.15. A probability can be read from Flgura 2-3 for any desired number of days 
in a sequence of favorable days after an unfavorable day. 


Figure 2-3 ia elso useful to an analysis of runs of unfavorable weather. 
If one assumes the present weal »er is good, for example, then the probability 
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of a change to inclement weather for one or more days is 0.14, approximately. 
Conditional probabilities of runs of unfavorable davs with greater minimum 
length are readily obtained from this graph. 

The final forecasting problem investigated was the prediction of ASTP 
launch conditions for a few days ahead, following a base day with known present 
weather and past weather. Further application of Markov theory in the form 
of the Chapman-Kolmogorov Equations was made, there being evidence of feasible 
predictions up to four days, at which time the '■•opirical value of p(U^|U Q ) 
generally reaches the unconditional value of p(U). These results are available 
in tabular form (Table 3-1) and also graphical form (Figures 3-1, 3-2). 

The theoretical results have been compared with the experimental counts 
to determine the amount of agreement between Markov chains of the order indicated 
by preliminary investigation of the first forecasting problem, and the data. 
Agreement was within 5 percent absolute value in all of the computations 
(Table 3-1). 

Finally, it is noted again that the definition of "unfavorable" used 
throughout this study is "the occurrence of one or more of the current ASTP 
constraints (subsection 2.1) at any hour between 1200 EST and 1700 EST". There 
is a known probability that an ASTP restriction will take place at a particu- 
lar hour on an unfavorable day in July and August, and as was demonstrated 
in subsection 2.5, this provides a link between the tabular values of this 
report and calculations of probabilities for specific hours. 
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